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APPLICATIONS  OF  THE  GENERALIZED  GLOBAL  EQUIVALENCE  RATIO  MODEL  (GGERM)  FOR 
PREDICTING  THE  GENERATION  RATE  AND  DISTRIBUTION  OF  PRODUCTS  OF  COMBUSTION 
IN  TWO-LAYER  FIRE  ENVIRONMENTS  - METHANE  AND  HEXANES 

Leonard  Y.  Cooper 

National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  02899 


ABSTRACT 

The  Generalized  Global  Equivalence  Ratio  Model  (GGERM)  was  developed  to 
predict  the  generation  rates  of  oxygen,  fuel,  and  other  products  of  combustion 
in  rooms  containing  fires.  The  GGERM  extends  to  general  transient  conditions 
the  global  equivalence  ratio  model  established  during  times  of  steady- state  in 
experimental  studies  involving  two-layer  compartment  fires.  The  present  work 
uses  the  GGERM  to  predict  the  time -dependent  upper- layer  mass  fractions  of 
products  of  combustion  (fuel,  oxygen,  CO,  and  others)  in  these  two -layer  fire 
experiments.  All  predicted  results  are  found  to  be  plausible  and,  where 
transient  data  are  available,  predicted  and  measured  results  compare  favorab- 
ly. However,  available  data  are  limited  and  additional  validation  of  the 
GGERM  under  more  varied  fire  conditions  is  required  before  it  can  be  used  with 
confidence  in  two-layer  zone-t)rpe  compartment  fire  models. 


Keywords:  building  fires;  combustion;  compartment  fires;  computer  models; 

equivalence  ratio;  fire  models;  global  equivalence  ratio; 
mathematical  models;  zone  models. 


INTRODUCTION  - THE  EXTENDED  UPPER  LAYER  IN  A ROOM  CONTAINING  A FIRE 


Background 

This  paper  presents  results  of  applying  and  validating  a general  model,  usable 
in  two- layer  zone -type  compartment  fire  analyses,  for  predicting  the  genera- 
tion rate  and  accumulation  of  combustion  products  throughout  a multi -room 
facility.  The  model  is  called  the  Generalized  Global  Equivalence  Ratio  Model 
(GGERM)  [1]^.  In  this  work,  the  GGERM  is  used  to  simulate  the  mass  fractions 
of  products  of  combustion  that  developed  during  the  transient  portions  of  the 
two -layer  fire  experiments  reported  in  references  [2]  and  [3] . Model 
predictions  are  presented  and  they  are  compared  to  experimental  results  when 
available . 

Consider  a multi-room  compartment  fire  and  assume  an  upper- layer/lower- layer 
zone -type  of  description  of  the  environment  in  each  of  the  rooms.  As 
indicated  in  Figure  1,  for  the  purpose  of  describing  combustion  processes  in 
any  room,  the  GGERM  defines  the  fire  and  plume,  which  may  originate  in  the 
lower  layer,  as  part  of  the  upper- layer  zone.  Also,  the  volumes  of  the  fire 
and  plume  are  assumed  to  be  negligible  compared  to  the  total  volume  of  the 
room.  Thus,  the  GGERM  treats  all  combustion  processes  in  any  room  with  a 
fire  as  occurring  in  the  upper  layer.  The  actual  combustion  process  is 
modeled  as  if  it  occurs  in  a well-stirred  reactor  with  a relatively  rapid 
throughflow  and  a correspondingly  negligible  residence  time. 


An  Overview  of  the  GGERM 

Presented  here  is  a brief  summary  of  the  essential  features  of  the  GGERM,  the 
details  of  which  are  presented  at  length  in  [1]. 

Let  Pj|.  represent  the  total  mass  of  a combustion  product  k in  the  upper 
layer  of  any  room  containing  a fire  during  an  arbitrary  multi -room  fire 
scenario.  The  goal  of  the  GGERM  is  to  estimate  the  combustion  contribution  to 
the  instantaneous  value  of  dPj^  u/^^»  rate  of  change  of  product  k in  the 

upper  layer. 

Conservation  of  species  leads  to  the  following  general  representation  for 

dPk.u/dt 

dp,  „/dt  - (1) 

where  the  right-hand- side  term,  Pj^  u»  represents  the  net  rate  of  product  of 
combustion  k flowing  to  the  upper  layer.  This  term  is  the  sum  of  material 
transfers  of  product  k to  the  extended  upper  layer  from  all  plumes,  jets. 


^Numbers  in  brackets  refer  to  the  list  of  REFERENCES  at  the  end  of  this 
report . 

list  of  NOMENCLATURE  is  included  at  the  end  of  this  report. 
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near -boundary  flows,  combustion  zones,  and  other  isolated  or  distributed 

sources . 

The  material  flowing  into  and  out  of  the  extended  upper  layer  generally 
includes  various  amounts  of  the  combustion  products  k associated  with  a fuel 
of  interest  and  other  inert  flow  components.  The  total  rates  of  mass  inflow 
and  outflow,  nipLowiN  %lowout  » corresponding  individual  rates  of 

inflow  and  outflow  of  products  k,  Pk,FLOWiN  ^k,FL0W0UT>  respectively,  are 

established  from  a variety  of  compartment  fire  modeling  considerations. 

These  rates  are  assumed  to  be  known. 

The  combustion  contributions  to  ^ are  modeled  by  the  GGERM.  As  depicted  in 
Figure  2,  the  GGERM  has  two  components:  the  accumulator,  which  simulates  spa- 
tially-averaged properties  of  the  generally  unsteady  upper  layer  environment; 
and  the  quasi -steady  reactor,  which  simulates  the  actual  combustion  processes 
taking  place  there. 

As  indicated  in  Figure  2 , all  of  the  flow  into  the  layer  is  modeled  as 
immediately  "entering”  the  reactor  where  it  may  participate  in  the  reaction 
process.  As  depicted  in  the  figure,  the  model  also  includes  a "feedback  loop" 
which  brings  flow  from  the  bulk  upper  layer  atmosphere  into  the  reactor  and 
the  reaction  process.  A critical  feature  of  the  overall  model  is  the  means  of 
evaluating  the  feedback  mass  flow  rate,  &peedbk • This  is  presented  in  Eqs . 
(42),  (43),  (46),  (47),  and  (49)  of  [1]. 

The  global  equivalence  ratio,  <i> , is  introduced  as  a normalized  ratio  of  the 
rates  of  fuel  to  oxygen  entering  the  reactor.  4>  is  defined  and  determined  by 


^ " ^^FUEL  , REACTIn/^02  , REACTIN  (2) 

where  PpuEL  reactin  ^02  reactin  mass  flow  rate  of  FUEL  and  02, 

respectively,  flowing  into  the  reactor;  Pj^  reactin  generally  defined  as  the 
mass  flow  rate  into  the  reactor  of  these  and  any  other  products  k;  and  r is 
the  known  fuel- to -oxygen  ratio  for  idealized  stoichiometric  combustion  of  a 
fuel  of  interest.  Note  that  the  ^ definition  of  Eq.  (2)  is  consistent  with 
the  4>  definitions  of  references  in  the  list  at  the  end  of  this  work. 

According  to  the  GGERM,  , the  rate  of  generation  of  an  arbitrary  product  k 
due  to  combustion,  is  given  by 


‘^/^FUEL  , REACTIN  “ 

where  the  fj^  functions  for  the  fuel  would  be  determined  from  steady  state 
experiments  of  the  type  carried  out  in  references  [2] -[9]. 

Having  obtained  values  for  , the  right-hand  side  of  Eq . (1)  can  be  evaluated 
from 


k . U 


k , FLOWIN 


k , FLOWOUT 


(4) 
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and  integration  of  Eq . (1)  can  proceed.  This  leads  to  solutions  for  the 


INITIAL  VALUE  PROBLEM  FOR  THE  MASS  FEIACTION  OF  FUEL,  OXYGEN,  AND  OTHER 
PRODUCTS  OF  COMBUSTION 

In  [1],  the  GGERM  implementation  of  Eq.  (1)  is  also  formulated  in  terms  of 
y , the  upper- layer  mass  fractions  of  products  k,  instead  of  Pj^  ^ . These 
are  defined  by 


, u ~ ^k  , u 


(5) 


where  m^  is  the  instantaneous  total  mass  of  material  in  the  upper  layer. 

The  governing  equations  for  Cp^EL  u > ^oz  u » other  Cj^  ^ of  interest  are 

[1]: 


^‘^^FUEL  , ^FUEL.U^FEEDBK^FUEL^*^)  " ^FUEL,U*FLOWIN 

^FUEL  , FLOWIN  ^FUEL^'^^^ 

^02  , U*FEEDBK^^02  ' ^02.U*FL0WIN 

^02,FL0WIn[^ 

^k,U*FL0WIN  (^FUEL  , U*FEEDBK  ^F  U E L , F L 0 W I N ^ ^k  ^ 

^^02  , U*FEEDBK  ^0  2 , F L 0 W I N ^ ^^k  ^ 

-=  (t)  (10) 

where  initial  conditions  for  CpypL  u»  ^02,  U’  ‘^k . u specified. 

Following  Figure  2 and  the  above  GGERM  description,  the  value  of  <f>  in  Eqs . 
(6)-(9)  is  calculated  from  Eq.  (2)  to  be 

^ t^^FUEL.FLOWIN  ^FUEL,U^FEEDBK^/^^0  2,FL0WIN  ^0  2 , U *^F  E E D B K ^ ^ 

(11) 

Also,  the  value  of  ®feedbk » calculated  from 

*FEEDBK  ^ fi^FEEDBK^l^y®’^  inflow  Conditions,  CpuEL.u*  ^02,0) 


niy  y ^ 


4 


which  represents  the  generalized  equation  set  presented  in  [1]. 

Finally,  the  functions  and  r are  assumed  to  be  known  and  specified 
properties  of  a fuel  of  interest. 

In  specific  fire  scenarios,  a closure  to  the  above  problems  for  the  Cj^  ^'s 
requires  specification  of  the  time -dependent  values  of  m^ , *plowin ’ 

PpuEL  FLOwiN  » ^02  FLOWiN  ’ Experimental  data  or  analytic  models  would  be 

expected  to  provide  these  latter  input  specifications. 

As  can  be  seen  from  Eqs . (6)  and  (7),  for  a particular  fire  scenario  the 
initial  value  problems  for  Cpyp^  u ^02  u coupled,  but  independent  of 

the  initial  value  problems  for  the  other  Cj^  u ■ solving  a given  problem, 

first  solutions  for  Cpup^  u ^02  u be  obtained.  Then  the  Aj^(t)  of 

Eq . (10)  would  be  determined  from  the  representations  of  the  right-hand  sides 
of  Eqs.  (8)  or  (9).  Finally,  the  Cj^  ^ for  other  combustion  products  of 
interest  would  be  obtained  from  solutions  of  Eqs.  (10). 

FIRE  SCENARIOS  WHERE  ALL  OXYGEN  INFLOW  IS  FROM  AN  AMBIENT  ATMOSPHERE 

Experiments  of  [2] -[10]  involved  test  configurations  depicted  in  Figures  3(a) 
and  3(b).  In  these  tests  all  components  of  Pq2  flowin  from  a standard 

atmosphere,  i.e.,  with  an  oxygen  mass  fraction  of  0.232.  Unless  noted 
otherwise,  analyses  and  results  throughout  the  remainder  of  this  work  are 
generally  applicable  only  for  fire  scenarios  where  Pq2  flowin  similarly 
from  a standard  atmosphere.  In  such  scenarios  it  can  be  shown  from  [1]  that 
Eqs.  (6) -(11)  can  be  written  as 

/^FLOWIN  ^^^FUEL  , “ ^F  U E L , U t ^F  U E L ^ ^*^F  E E D B K /*F  L 0 W I N ^ 

-H  0.232r<^(ss)  [1  + ^^  (^)  ] / [ l + 0 . 232r(?!.(  ss  ) ] 

(13) 


(mu/ftlpLOWIN  ^^*^02  , = Cq2  ^ u [^*^^02  ^’^FEEDBk/*FLOWIN  ) ■ 

0.232[1  + + 0.232r<^(SS)  ] 

(14) 


(%/*FL0WIN  + Cj^  u “ t ^FUEL  , U ^*FEEDBk/%L0WI  N ^ 

-h  0.232(ji(SS)r/[l  -t-  0.232r<?i.(ss)  ] )f^  (^) 

(15) 


“ { ^02 , u ^*feedbk/*flowin ^ 

-H  0.232/[l  -H  0.232r.^(ss)  ] (,ji)  (16) 

»A,(t)  (17) 
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^ + (’^FEEDBk/^FLOWIN  ^‘^FUEL  . U t 

{1  + (ftpEEDBK/^'FLOWI  N )^02  , U 


1 + l/(0.232r(?i(ss))])/ 

^(SS)j.[i  + l/(0.232r<^(ss)  )]  ) 

(18) 


As  in  [1],  in  Eqs . (13)-(18)  is  defined  as  the  value  of  ^ under  virtual 

steady-state  conditions,  i.e.,  the  value  of  </>  when  ^feedbk  ""  ^ when 

Pfuel  flowin  ^02  FLOWiN  maintained  at  their  current  values,  i.e.. 


^(SS)  = (p 


FUEL.FLOWI 


n/^0 


2 . FLOWI 


n)A 


(19) 


As  is  also  defined  in  [1],  (<^^  ^^  ^ ) is  the  product-k  mass  fraction 

flowing  out  of  a GGERM  reactor  under  steady  state  conditions  when  the  Global 
Equivalence  Ratio  of  the  inflow  to  the  reactor  is  . 

As  presented  in  [1] , for  any  product  k there  is  a one-to-one  correspondence 
between  f ^ actoui 

clgitREACioui  = 0.232^(ss,r[l  + (^(  ss ) ) ]/(i  + 0.232^<ss>r)  (20) 


cil!;REACioui  - 0.232[1  + r,^<ss)f^^(^(ss))]/(i  + o.  232r,^<  = = > ) 


(21) 


c(  ss  ) 

'-k  , REACTOUT 


0 . 232r^^  ^ fj^  (0^  0/(1  + 0.232r<^^2^O  fot  k other  than 

fuel  or  oxygen 

(22) 


Solving  for  the  fj^  in  the  above  leads  to 


^FUEL  [1  1/(0  • 232r^^  ss  ) ) j c|,ss^  ^ ss  ) ) . p (23) 

^02  “ ^oi^^EACTOUT 

[c^i"^EACTOUT(^^""M  - 0.232]/(0.232r^(ss))  (24) 

fj^  - [1  + l/(0.232r(^(ss))]c^sy^^^^y^(<^(ss))  for  k other  than  (25) 

fuel  or  oxygen 

When  data  for  the  (<^^  ^^  O functions  for  a particular  fuel  are 

presented  in  reports  of  Figure-3- type  experiments  in  a standard  atmosphere,  it 
is  possible  to  determine  analytic  representations  of  these  and  to  then 
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construct  the  fj^  functions,  required  for  use  in  the  GGERM,  by  using  Eqs . (23)- 
(25).  Similarly,  known  fj^  functions  can  be  used  to  construct  the 
by  using  Eqs.  (20) -(22). 

Eqs.  (13) -(15)  are  completed  with  specification  of  ftipEEDEK/’^FLOWiN  • fire 

scenarios  considered  here,  this  is  found  from  [1]  to  be 


*FEEDBk/*FLOWIN  ™^^^*FEEDBK  , 1 /*FL0WI  N ’ ^FEEDBK  , 3 /*^FLOWI  N ^ (26) 


where 


fti. 


F E E D B K , 1 /*F  L 0 W I N 


0 if  1 or  ^02. U ■ ^02?REACT0UT  “ 

[(1  - <?^fuel/^'""')/(<?^fuel/^'""'  + 0-232r<?ipuEL)]- 

[(”0  - c(SS)  W('0  232  - r(SS) 

‘■'^0  2,U  ^0  2 , RE  AC  T 0 U T ^ ''0  2,REACT0UT'-I 

if  ^^^^V^fUEL  ^ ^ ^0  2 U ■ ^oi^REACTOUT  ^ ^ 

(27) 


fil 


FEEDBK , 3 /*FLOWIN 


^ 0 if  < 1 or  CpuEL.U  - ^FUEL.REACTOUT  “ 

- l)(0.232r<^(ss))/(i  + o.232r^(ss))]  . 

t^^fUEL.U  " ^^UeZ  , REACTOUT  ‘ U E Z , RE  A C T 0 U T ^ ^ 

rf  ^ ^ ^FUEL.U  ■ uiZ  , REACTOUT  ^ ® 

(28) 


and  where,  gs^  ^ and  c^I^^eactout  obtained  from  Eqs.  (20)  and 

(21),  respectively. 

In  the  above,  </>puel  defined  as  in  [1]  as  the  maximum  value  of  4>  which  leads 
to  complete  combustion  of  the  fuel,  i.e.,  the  maximum  value  of  4>  where 

^FUEL  ■ ^FUEL, REACTIN'  Using  Eq . (3) 

<^fuel  ^ maximum  root  of  fFUEL  ^*^fuel  ^ -f  1 = 0 (29) 


It  is  assumed  that  <^fuel  ^ ® bounded. 

Similarly,  4>q2  is  defined  as  in  [1]  as  the  minimum  value  of  which  leads  to 
complete  consumption  of  oxygen,  i.e.,  the  minimum  value  of  4>  which  leads  to 
i*0  2,REACTiN  “ Using  Eq.  (3) 
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<^Q2  = minimum  root  of  4>o2^0  2^^0  2^^  + 1 = 0 


(30) 


It  is  assumed  that  <^q2  >0  and  is  bounded. 

CONSTRUCTING  ANALYTIC  REPRESENTATIONS  FOR  THE  PRODUCT  GENERATION  FUNCTIONS  AND 
THE  STEADY  STATE  MASS  FRACTIONS  FROM  EXPERIMENTAL  DATA 

General  Considerations 


In  general,  application  of  the  GGERM  for  a given  fuel  requires  calculations 
based  on  analytic  representations  of  its  f ' s and  ^ ' s . These  would 

be  obtained  from  steady  state  experiments  of  the  kind  described  in  references 
[2],  and  [4]-[9]. 

If,  as  in  [2],  [4],  [5],  and  [9],  reduced  data  points  are  provided  for  the 
^k^REACTOUT ^ ® » then  these  data  points  must  be  used  to  construct  analytic, 
curve -fitted  representations  of  these  functions.  These  latter  representations 
would  be  used,  in  turn,  to  calculate  a consistent  set  of  the  f ^ ' s . The 
calculation  would  use  Eqs . (23) -(25). 

If,  as  in  [7]  and  [8],  reduced  data  points  are  provided  for  the  f j.  , then  these 
data  points  can  be  used  to  construct  analytic  curve -fitted  representations  of 
these  functions.  These  latter  representations  would  be  used,  in  turn,  to 
calculate  a consistent  set  of  the  ^ ® this  case,  the  calcula- 

tion would  use  Eqs.  (20) -(22). 

The  curve-f itting  task  must  be  carried  out  so  that  the  resulting  fj.  (^)  and 
^k ^REACTOUT ^ ^ functions  are  physically  and  analytically  well-behaved  for 
arbitrary  non-negative  values  of  (f>  and  . The  following  mutually- 

consistent  physical  and  analytic  constraints  are  required: 

1.  The  fj^(^)'s  and  ^ (^^  ss  ) ) » g gj-e  continuous  functions  for 

arbitrary  non-negative  values  of  <j>  and 

2.  According  to  the  Eq.-(29)  definition  of  ^fuel  , 0 < ^ < ^fuel 

0 in  particular,  all  fuel  entering  the  reactor  is  completely 
consumed,  i . e . , 


LO 


FUEL 


- P 


FUEL , REACT  IN 


for  0 < <f>  :£  4> 


FUEL 


(31) 


Consistent  with  whatever  reactor  chemistry  is  applicable  as  </>-*•  0,  the 
generation/consumption  rate  of  oxygen  and  all  products  k will  be 
directly  proportional  to  the  generation/consumption  rate  of  fuel.  The 
constants  of  proportionality  are  designated  in  terms  of  the  characteris- 
tic product-k  constants,  q » 
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lim  0\  = - 
4>  0 


1 


, 0 U E L 


(32) 


FUEL  , 0 


(33) 


where 


^02  0 complete  stoichiometric  combustion  at  (/»  -=  0 

(34) 

Using  Eq.  (3)  and  (32),  with  k -♦  O2  (i.e.,  k representing  O2 ) , and  Eq. 
(31),  and  then  using  Eq . (21)  leads  to 


lim  „ (35) 

<!>  -*  0 


iJ.ui  '-OXY  , REACTOUT  > 

<^(SS)  0 

0.232[1  + . 0. 232^^2, q)  + 0(<jicss)2)] 

(36) 

Using  Eq.  (3)  with  k ->  FUEL  and  Eq.  (31),  and  then  using  Eq.  (20)  leads 
to 


^FUEL^^^  ^ ® — ^ — ^fuel  (37) 


c(ss)  /^(SS)\  « 0 for  0 < m ('38) 

Using  Eqs . (3)  and  (32)  with  k other  than  FUEL  or  O2 , and  Eq.  (31),  and 
then  using  Eq.  (22)  leads  to 


For  k other  than  FUEL  or  O2 : 

llm  (^)  - 0 (39) 

<f>  -*  0 


lim  ^ “ 0 . 232r<^^  ^ Fj^  Q 

<^(ss)  _ 0 


(40) 
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3. 


According  to  the  Eq . (30)  definition  of  ^02 » ^ — ^02  ^ 

(^  ->•  00  in  particular,  all  oxygen  entering  the  reactor  is  completely 
consumed,  i . e . , 


02 


= - P 


02 , REACTIN 


when  (f>  > 4>q  2 ^ ® 


(41) 


Consistent  with  whatever  reactor  chemistry  is  applicable  as  ^->00,  the 
generation/consumption  rate  of  oxygen  and  all  products  k will  be 
directly  proportional  to  the  generation/consumption  rate  of  fuel.  The 
constants  of  proportionality  are  designated  in  terms  of  the  characteris- 
tic product-k  constants,  ,«,  • 


lim  Coy. 

(f>  -*  CO 


. CO  U E L 


(42) 


r 


FUEL  , CO 


= i 


(43) 


where 


^02  CO  “ ■ complete  stoichiometric  combustion  at  0 -*•  <» 

(44) 

Using  Eq . (3),  with  k O2 , and  Eqs . (2)  and  (41),  and  then  using  Eq. 
(21)  leads  to 


^02^*^)  “ ■ V(i^<^)  when  > 4>q^  > 0 (45) 

n(.  ss  ) ( ^ ^ ^ ^ 0 when  ^ ^ ^ ^ >■  <i>n « ^ 0 (46) 

Using  Eqs.  (3)  and  (42),  with  k -*  FUEL,  and  Eqs.  (2)  and  (41),  and  then 
using  Eq.  (20)  leads  to 

lim  ffuEL^*^)  = when  <f>  > > 0 (47) 

(p  ao 

lim  “ 

^ ( S S ) -4.  00 

1 + (1  - ro2  ./0.232)[l/(r^<ss>ro2  .)]  + 0(l/.#css)2) 

(48) 

Using  Eqs.  (3)  and  (42)  with  k other  than  FUEL  or  O2 ^ Eq.  (2)  and  Eq. 
(41),  and  then  using  Eq.  (22)  leads  to 
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(49) 


For  k other  than  FUEL  or  O2 : 


lim  = - rj^,«/(r<^ro2,a.) 


1 i_ni  p ( s s ) 

J.XU1  '-PUEL.REACTOUT 
0 ( S S ) 00 


(0( SS ) ) 


r^,./(r^!ss)r„2  „) 


(50) 


As  discussed  in  [1],  when  reaction  processes  only  involve  complete  stoichiome- 
tric combustion,  </'fuel  » ^02  ’ given  by 


^(STOICH)  = M(STOICH) 
‘^FUEL  ^0  2 


= 1 


f(ST0ICH)(^) 


' if  0 < (^  < 1 

^ rj^/<;6  if  <f>  > 1 


where 


FUEL 


1.  r 


02 


l/r 


(51) 


(52) 


(53) 


and  where  the  other  Fj^  can  be  determined  from  the  stoichiometric  reaction  of 
the  particular  fuel  of  interest.  These  results  are  consistent  with  the 
results  of  items  1-3  above. 


The  Product  Generation  Functions  and  the  Steady  State  Mass  Fractions  for  the 

Combustion  of  Methane  Using  Data  of  [21 

Reference  [2]  provides  reduced  steady  state  combustion  of 

natural  gas  in  the  test  configuration  of  Figure  3b.  These  data  were  acquired 
over  a wide  range  of  layer  temperatures,  but  these  never  exceeded  800K. 
values  were  in  the  range  0 < <3,  in  these  tests,  in  addition  to  the 

rate  of  oxygen  supply  to  the  extended  upper  layer  through  entrainment  into  the 
burner's  fire  plume,  the  gas  in  the  upper  layer  collector  was  injected  with  a 
specified  flow  of  air  and  its  associated  oxygen.  Some  of  the  reference- [ 2 ] 
test  results  and  results  of  other  natural  gas  tests  using  earlier  versions  of 
the  same  test  configuration  were  presented  in  [4],  [5],  and  [9]. 

It  will  be  assumed  here  that  the  natural  gas  fuel  used  in  the  experiments  of 
[2]  can  be  modeled  as  pure  methane,  CH^ . 

c^ss^actout  ^ ' ^2’  ^^2  » ^2 ^ ^2  presented  in  [2] 

and  are  reproduced  here  in  Figures  4-10.  Analytic  approximations  for  the  fj, 
and  functions  which  characterize  the  combustion  of  CH^  will  be 

constructed  by  using  these  data,  the  analytic  constraints  of  the  above  items 
1-3,  and  Eqs.  (23) -(25).  With  these  functions,  the  GGERM  will  be  used  to 
simulate  the  transient  upper- layer  environment  in  the  methane  combustion 
experiments  of  [2]. 
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It  is  important  to  point  out  that  the  data  of  the  reference -[ 5 ] experiments  of 
Toner,  included  in  Figures  4-9,  were  acquired  at  higher  layer  temperatures, 
greater  than  900K,  than  were  all  the  other  data  in  the  figures.  Regarding  the 
generation  of  CH^  , and  especially  O2 , CO2 , and  H2O,  Figures  4-7  show  that  at 
these  temperatures  the  reactions  closely  simulate  complete  stoichiometric 
combustion.  The  Toner  data  suggest  that  GGERM  simulations  at  higher  tempera- 
tures would  be  improved  if  the  f ' s of  Eq . (3)  were  generally  assumed  to  be 

dependent  on  layer  temperature  as  well  as  <f) . For  now  this  temperature 
dependence  is  ignored  and  IN  ALL  CASES  THE  DATA  POINTS  OF  THE  REFERENCE -[ 5 ] 
EXPERIMENTS  OF  TONER,  INCLUDED  IN  FIGURES  5-10,  WILL  BE  DISREGARDED  IN  THE 
ANALYSIS  TO  FOLLOW.  Further  discussion  of  this  will  be  presented  toward  the 
end  of  this  work. 


Methane  (CH^ ) : 

The  reduced  data  for  ^ is  presented  in  Figure  4 (a  re-labelled 

version  of  Figure  4.5  of  [2]).  Included  in  the  figure  is  a plot  of 
^CH 4 ^ REACTOUT  Complete  stoichiometric  combustion.  The  latter  plot  is 

obtained  from  Eq . (20)  using  Eqs . (51) -(53)  and  r = 1/4. 

From  visual  observation  of  the  data,  the  value  of  </'puel  *^ch4  approximated 
as 


^FUEL  “ <^CH4  = 


(54) 


It  is  now  observed  from  Figure  4 and  from  the  results  of  the  items  1-3  that 
for  CH^ , and  possibly  for  other  fuels,  certain  basic  functional  characteris- 
tics of  the  REACTOUT  approximately  similar  for  real  combustion  and 

for  complete  stoichiometric  combustion,  viz.,  reactout  identically 

zero  for  0 < ^/•^puel  - monotonically  increasing  with  uel 

<^/(^fuel  ^ increasing  with  a positive  slope  (designated  here  as  a)  from  zero 
at  '^/5'puel  as3niiptotic  to  1 as  ^/<^puel  "*■ 

All  of  the  above  characteristics  are  shared  by  the  following  two-parameter 
representation  which  satisfies  the  required  asymptotic  behavior  of  Eq.  (48) 


c 


( ss ) 

FUEL , REACTOU 


T 


= [1  - 


(^FUEl/^^  ^ ^ [1/(^*^FUEL^  “ ^ n<^FUEL ^ 


if  <?^/<^FUEL  > f (55) 


where 


a = lim  feactout/^*^  (56) 

•^/^FUEL  i i 

Eq.  (55)  is  exact  in  the  case  of  complete  stoichiometric  combustion  where, 
according  to  Eqs.  (20)  and  (47) -(49),  <^puel  ^ ^ ^ i/[i  1/(0. 232r)]. 

Indeed,  it  is  precisely  the  j^p^CTOUi  °f  • (55),  with  the  latter  values 

of  <^puel  ^^*5  o,  that  is  plotted  in  Figure  4. 
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data . 


The  Eq.  (55)  representation  is  used  here  to  fit  the  ^ 

From  the  Figure  4 data,  o is  approximated  as 


a - 0.037 


(57) 


The  desired  curve -fit  is  therefore  taken  to  be 


c(  ss  ) 

, REACTOUT 


0 if  </)  < 0.70  <^fuel  *^cha 

[1  - (0.70/<^)  ]/{l  + [1/(0.037-0.70)  - l](0.70/<^)) 

if  </>  > 0.70  = <^PUEL  “ ^CHA 


(58) 


This  is  plotted  in  Figure  4 


The  reduced  data  for  c^I^reactout  presented  in  Figure  5 (a  re-labelled 
version  of  Figure  4.3  of  [2]).  Included  is  a plot  of  for 

complete  stoichiometric  combustion.  This  is  obtained  from  Eq.  (21)  using  Eqs . 
(51) -(53)  and  r = 1/4 


c(  ss  ) 

*-02  . REACTOUT 


STOICH 


0.232(1  - <t>)/{l  + 0.232^/4)  if  0 < ^ < 1 
0 if  1 < 0 


(59) 


As  can  be  seen  in  the  figure,  when  0 < <^  < 0.60  an  approximation  of  complete 
stoichiometric  combustion  provides  a good  estimate  for  the  real  combustion 
data.  This  approximation  will  be  used  for  the  desired  curve -fit  in  this  range 
of  4> . 

From  extrapolation  of  data  in  Figure  5,  4>q2  approximated  as 


4>q2  = 3.50  (60) 

In  the  range  0.60  < 4>  < 3.50  the  function  A (3. 50  - with  A chosen  to 

force  continuity  with  Eq . (59)  at  (/>  = 0 . 6 , is  found  to  fit  well  the  Figure  5 
data . 

Consistent  with  the  above,  the  desired  jqut  curve-fit  is  taken  to  be 
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c(  SS  ) 

'“0  2,REACT0UT,ST0ICH 


if  0 < (^  < 0.60 


{ (0.232) (0.40)/[l  + (0.232) (0.60)/4] ) • 

[(3.50  - (f>)/2.90]/{4>/0.^0)^'^ 

if  0.60  < <f)  < 3.5 

^02  REACTOUT  STOJCH  ~ ^ 3.50  ^ (f) 

(61) 


This  is  plotted  in  Figure  5. 


CO2  and  H2  0 : 

The  reduced  data  for  and  ^ are  presented  in  Figures 

6 and  7 (re-labelled  versions  of  Figures  4.7  and  4.8  of  [2]),  respectively. 
Included  are  plots  of  these  for  complete  stoichiometric 

combustion.  These  are  designated  as  c^ffEACTOUT  stoich  they  are  computed 

from 


c(  SS  ) 

'"C  02, REACTOUT,  STOICH 


(0.232/4)  (ll/4)<^/(l  + 0.232<^/4)  if  0 < 0 < 1 
(0.232/4)(ll/4)/(l  -I-  0.232<^/4)  If  1 < 4> 


(62) 


c(  SS  ) 

'-'H  20, REACTOUT,  STOICH 


(9/ll)c(SS) 


REACTOUT , STOICH 


(63) 


Eqs . (62)  and  (63)  were  obtained  from  Eq.  (22)  using  Eqs . (52),  r = 1/4,  and 
Fco2  = 11/^  and  rH2o  “ 9/4  [1]. 

As  can  be  seen  in  Figure  6,  when  0 < <^  < 0.60,  a complete  stoichiometric 
combustion  approximation  provides  a good  estimate  for  the  real  combustion  data 
of  Ccoi ^ REACTOUT  • simplicity,  this  approximation  will  be  used  here  for 

the  somewhat  extended  range  0 < 4)  < 0.70.  Then,  for  the  range 
0.70  < 4>  ^ 3.00,  ^ will  be  approximated  as  being  constant. 

Finally,  consistent  with  the  required  asymptotic  behavior  of  Eq . (49), 

^coi  ^ REACTOUT  approximated  as  being  proportional  to  l/<t>  for  4>  > 3.00. 

Consistent  with  the  above,  and  with  the  required  continuity  of  the 
approximating  function,  the  desired  ^ curve-fit  is  taken  to  be 
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r ( SS  ) 

'-C02,REACT0UT 


0 2 ^ R E A C T 0 U T , S T 0 I C H ^ ~ ^ ~ 

" Aco2  1 if  0.70  < <i>  < 3.00 
(3.00/(?!))Aco2,  1 if  <^  > 3.00 


(64) 


where 


^co2.i  = (0.232/4)(ll/4)(0.70)/[l  + (0 . 232) (0 . 70)/4]  = 0.107...  (65) 

This  is  plotted  in  Figure  6. 

As  can  be  seen  in  Figure  7,  when  0 < <^  < 0.70  and  2.25  < <^  it  is  reasonable  to 
use  the  complete  stoichiometric  combustion  approximation  to  estimate  the  real 
combustion  data  of  ^ . In  the  range  0.70  < 4>  < 2.25  the  fitted 

function  is  completed  by  the  straight-line  segment,  plotted  in  the  figure, 
that  joins  the  curves  applicable  in  the  two  outer-(^  ranges. 

Consistent  with  the  above,  the  desired  ^ curve-fit  is  taken  to  be 


r(  SS  ) 

*^820  , RE  ACTOUT 


c(.ss)  if  0 < < 0 70 

^H20,REACT0UT,ST0ICH  ~ ^ 

^ ^E20.1<?^  " ^H20,2  if  0-70  < < 2.25 

c(  SS  ) if  2 25  < (i 

''H20,REACT0UT,ST0ICH  ^ ^ ^ 


where 

^B20,1  - - c^iE;i)/(2.25  - 0.70) 

A,2o,2  “ (0-70c<|e;,  - 2.25c^|g>i)/(2.25  - 0.70) 
=H2o!i  “ (0.232/4)(9/4)(0.70)/[l  + 0 . 232 (0 . 70)/4] 
c<|S)2  - (0.232/4)(9/4)/[l  + 0 . 232 (2 . 25)/4] 


(66) 


(67) 


CO,  and  H2 : 

The  reduced  data  for  and  Ch2?reactout  presented  in  Figures  8 

and  9 (re-labelled  versions  of  Figures  4.7  and  4.2  of  [2]),  respectively. 

As  seen  in  Figure  8,  the  data  show  that  is  effectively  zero  for  4> 

less  than  approximately  0.55.  For  4>  > 0.55  and  as  shown  in  the  figure,  the 
data  curve-fit  will  be  accomplished  with  straight-line  segments  in  the  ranges 
0.55  < 4>  < 1.15  and  1.15  < 4>  < 3.00.  Also,  consistent  with  the  required 
asymptotic  behavior  of  Eq . (49),  will  be  approximated  as  being 

proportional  to  l/<^  for  4>  > 3.00. 
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Consistent  with  the  above,  the 


curve -fit  is  taken  to  be 


0 if  0 < 0 < 0 . 55 


c(  ss  ) 

^CO,REACTOUT 


CO  , 1 

^CO  , 2 

if 

0.55 

< 4>  < 

1.15 

CO , 3^  ~ 

0 A 

if 

1.15 

A 

-0- 

lA 

3.00 

0. 0205(3. 00)/<ji  if  3.00  < 4> 


(68) 


where 

Aco, 1 = 0.0127/(1.15  - 0.55) 

\o,2  = 0.55(0.0127)/(1.15  - 0.55) 

(69) 

^co,3  = (0.0205  - 0.0127)7(3.00  - 1.15) 

A^o  ^ = [1.15(0.0205)  - 3.00(0.0127)]/(3.00  - 1.15) 


The  curve-fit  of  the  c^ffREACTOUT  of  Figure  9 is  done  with  three  curve 

segments,  in  a manner  analogous  to  that  of  c^o?reactout • resulting 

function,  plotted  in  the  figure,  is  taken  to  be 


REACTOUT 


0 if  0 < ^ < 0.90 

[0.00145/(2.00  - 0.90)]-^  - 0 . 90(0 . 00145)/(2 . 00  - 0.90) 

if  0.90  < <;6  < 2.00 

0.00145  if  2.00  < <f>  < 3.00  (70) 

0.00145(3.00)/^  if  3.00  < 4> 


The  Product  Generation  Functions  and  the  Steady  State  Mass  Fractions  for  the 

Combustion  of  Hexane  Using  Data  of  [61 

Reference  [6]  provides  reduced  data  for  steady  state  combustion  of 

hexane,  CgH^^,  in  the  test  configuration  of  Figure  3a. 

The  plotted  fj^  data  for  k , O2  , H2O,  CO2  , CO,  and  H2  presented  in  [6] 

and  are  reproduced  here  in  Figures  10-15.  Analytic  approximations  for  the 
fj^  (^)  and  ^ ) functions  which  characterize  the  combustion  of 

the  fuel  will  be  obtained  by  first  applying  the  analytic  constraints  of  the 
above  items  1-3  to  these  data  and  then  using  Eqs . (20) -(22).  With  these 
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functions,  the  GGERM  will  then  be  used  to  simulate  the  transient  upper- layer 
environment  in  the  hexane  combustion  experiments  of  [3]. 


for  k Cg  ^ , O2  , H2  0 , and  CO2  : 

The  reduced  data  for  f , k -*•  FUEL  (i.e.,  O2  , H2O,  and  CO2  , are 

presented  in  Figures  10-13  (re-labelled  versions  of  Figures  4-35,  4-34,  4-36, 
and  4-32  of  [6]),  respectively.  These  are  all  the  components  that  would  be 
involved  in  the  complete  stoichiometric  reaction  of  the  fuel.  Included  in 
each  of  the  figures  is  a plot  of  the  corresponding  fj^  function  for  the 
complete  stoichiometric  combustion  reaction.  The  latter  plots  were  obtained 
using  Eqs . (51) -(53).  Note  that 


For  : 

r - 43/152,  . . 1,  ro2  = - 152/43, 

Thzo  - 63/43,  r<.02  - 132/43 


(71) 


It  is  now  observed  from  Figures  10-13  that  for  CgH^^ , and  possibly  for  other 
fuels,  there  are  positive  values  of  <f> , denoted  by  stoichi  ^ such  that  in 
the  range  0 < <i>  < <f>^  stoichi*  bhe  present  f^.  can  be  well-fitted  by 

f^STOiCH)  of  Eq.  (52^.  Such  representations  would  satisfy  the  required 
criteria  of  Eqs . (35),  (37),  and  (39).  (Note  that  *^puel’Stoichi  ~ 
is  also  possible  to  choose  positive  values  of  (f> , denoted  by  <f>^  stoichz  * 
possibly  out  of  the  range  of  the  fj.  data  or  out  of  the  range  of  ^'s  expected 
in  real  two -layer  fire  scenarios,  such  that  in  the  range  4>  > 4>^  stoich2 
f^STOiCH)  of  Eq.  (52)  again  provides  representations  for  the  f^.  which  are  good 
fits  to  the  data,  if  they  exist  in  the  range,  or  are  at  least  consistent  with 
the  required  criteria  of  Eqs.  (45),  (47),  and  (49).  (Note  that  <^oxy »stoichi  = 
*02-') 


^6  A ’ ® » SUd  CO2  I 

Using  the  above  ideas,  choices  for  the  4>^  stoichi  stoich2  made  and 

the  following  analytic  representations  are  used  to  construct  (the  four)  fj^ 
functions  which  are  well -fitted  to  the  data  and  continuous  for  all  4> , and 
which  satisfy  the  required  criteria  of  Eqs.  (35),  (37),  (39),  (45),  (47),  and 
(49): 


^k  i f ^ ^ .STOICHI  ’ 

(Ik/^^t^  " " *^k^^‘^k,ST0ICH2  ” ^^/^^k,STOICH2  “ ^k.STOICHl^^ 

^k,  STOICHI  ^ ^ ‘^k.ST0ICH2’ 

< Ik/^  ^ - *^k,ST0ICH2 

(72) 
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The  following  values  of  <^k.STOicHi  ^k,ST0iCH2  chosen  for  the  four 

curve  fits : 


^C6H1  A , STOICHl  ^ ‘^FUEL  ^ 0-85 
^0  2,ST0ICH1  ""  0.80 
*^H20,ST0ICH1  ""  0.85 
*^002  .STOICHl  0.70 
‘^k,ST0icH2  ^ ^0  2 ^ 10.0  for  k 


(73) 


Cg  Hi  ^ , O2  , H2  0 , and  CO2 


Note  that  Eq . (72)  is  exact  for  complete  stoichiometric  combustion  in  which 

case  all  stoichi  ^ ^k  stoich2  ""  middle  range  of  <}>  shrinks  to  a 

point. 

The  fj^  representations  of  Eqs . (72)  and  (73)  are  plotted  in  Figures  10-13. 


CO , and  H2 : 

The  reduced  data  for  f^^Q  and  fg2  are  presented  in  Figures  14  and  15  (re- 
labelled Figures  4-33  and  4-37  of  [6]),  respectively. 

As  seen  in- Figure  14,  the  data  show  that  f^Q  is  effectively  zero  for  4>  less 
than  approximately  0.55.  For  ^ > 0.55  and  as  shown  in  the  figure,  the  data 
curve -fit  will  be  accomplished  with  straight-line  segments  in  the  ranges  0.55 
< 4>  < 1.00  and  1.00  < <f>  < 3.50.  Also,  consistent  with  the  required  asymptotic 
behavior  of  Eq.  (49),  will  be  approximated  as  being  proportional  to  l/<f> 

for  4>  > 3.50. 

Consistent  with  the  above,  the  fj,Q  curve-fit  is  taken  to  be 


0 if  0 < ^ < 0.55 


(43/84) fco 


0.200<^  - 0.110  if  0.55  < <^  < 1.00 
0.090  if  1.00  < <^  < 3.50 


0. 090(3. 50)/<^  if  3.50  < <t> 


(74) 


The  curve-fit  of  the  f^2  data  of  Figure  15  is  done  with  three  curve  segments, 
in  a manner  analogous  to  that  of  f^o • resulting  function,  plotted  in  the 
figure,  is  taken  to  be 
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0 if  0 < < 0.70 


(43/7)fH2 


0.050(<^/0.70  - 1)  if  0.70  < <t>  < 1.40 
0.050  if  1.40  < <f)  < 3.50 


0.050(3. 50)/<^  if  3.50  < 4> 


(75) 


Plots  of  the  Product  Generation  Functions  and  the  Steady  State  Mass  Fractions 

for  the  Combustion  of  Methane  and  Hexanes 


The  fj^  functions  for  CH^  were  computed  from  Eqs . (23) -(25)  using  the 

representations  of  Eqs.  (54) -(70).  The  f pair  for  each  of 
the  six  products  are  presented  in  Figures  16-21. 

The  functions  for  CgH^^  were  computed  from  Eqs.  (20) -(22)  using 

the  fj^  representations  of  Eqs.  (71) -(75).  The  f pair  for  each 

of  the  six  products  are  presented  in  Figures  22-27. 

All  the  above  fj.  , ^k^lEACTOUT  representations  will  be  useful  in  GGERM- 
implemented  two -layer  zone -type  compartment  fire  model  simulations  of  fire 
scenarios  involving  the  combustion  of  CH^  or  CgHj^^.  Results  of  such  simula- 
tions will  be  presented  in  the  next  section. 


SIMULATING  EXPERIMENTS  OF  REFERENCES  [2]  AND  [3] 

The  Combustion  of  Methane  in  Figure-3b  Test  Confignirations 

Consider  the  64  CH^  tests  reported  in  Appendix  B of  [2] . The  test  configura- 
tion for  these  is  depicted  in  Figure  3b.  Except  for  initial  test  time 
intervals  (on  the  order  of  1 to  2 minutes) , which  were  relatively  short 
compared  to  the  times  to  reach  steady  state  conditions  (on  the  order  of  20  to 
30  minutes),  the  layer  interface  elevation  in  all  of  these  tests  was  constant, 
slightly  below  the  bottom  of  the  upper  layer  collector.  In  the  simulations  to 
follow  it  is  assumed  that  the  elevation  of  the  layer  interface,  i.e.,  the 
volume  of  the  upper  layer,  is  constant  throughout  a test  run.  From  [2]  the 
volume  is  estimated  to  be  the  volume  of  the  collector,  3.89  m^ . 

Each  test  was  carried  out  with  a fixed- elevation  gas  burner  surface  and  with 
constant  upper- layer  supply  rates  of  fuel  and  injected  upper- layer  air. 

Because  of  the  constant- interface-elevation  assumption,  which  implies  a 
constant  rate  of  air  entrainment  from  the  lower  layer  (i.e.,  the  laboratory 
environment)  to  the  fire  plume,  each  test  is  simulated  with  constant  values  of 

*FL0WIN  » ^FUEL  . FLOWIN  > ^02,FL0WIN 

^FLOWIN  ^FUEL.FLOWIN  ^AIR.FLOWIN  ^FUEL.FLOWIN  ^0  2,FL0WIn/®'^^^ 

(76) 
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In  Eq.  (76),  Pair.flowin  defined  as  the  mass  rate  of  flow  of  air  into  the 
upper  layer. 

With  the  above  assumptions,  solutions  to  the  GGERM  Eqs . (13) -(30)  have  been 
carried  out  to  simulate  experiments  of  reference  [2].  For  each  experiment 
simulated,  calculations  were  carried  out  with  the  two  different  models  of  CH^ 
combustion,  viz.,  real  combustion  and  complete  stoichiometric  combustion  In 
the  simulations  it  was  assumed  that  the  initial  mass  fractions  of  products  of 
combustion  in  the  upper  layer  were  the  same  as  those  of  an  ambient  atmosphere, 
which  were  approximated  by  Cq2  ^ 0.232,  and  all  other  Cj^  ^(t  = 0)  = 

0. 

Results  from  Appendix  B of  [2]  are  presented  here  for  the  simulation  of  test 
run  1.  For  this  test  the  initial  and  characteristic  final  temperature  of  the 
layer  were  reported  to  be  init  ^ final  ^ 529.7  K,  respective- 

ly. In  the  simulation,  the  upper  layer  temperature  was  assumed  to  be  at  the 
constant  reasonable  value  (2/3)  (T^  pj^AL  " init^  451.8  K. 

The  mass  of  the  upper  layer  was  computed  according  to  a perfect  gas  law  model 
of  the  upper  layer  with  the  gas  constant  taken  to  be  that  of  a standard 
atmosphere.  The  other  conditions  of  the  test  were  taken  from  [2]  to  be: 


^FUEL.FLOWiN  0.00135  kg/s 
^AIR.FLOWIN  0.01030  kg/s 


(77) 


In  this  test  run,  all  inflowing  air  was  from  plume  entrainment.  Using  Eq. 

(77)  in  Eq.  (19),  the  constant  simulated  value  of  foj-  this  experiment  is 

found  to  be  2.33.  The  reported  in  [2]  was  2.17.  The  difference  in 

these  two  values  is  attributed  to  the  fact  that  in  the  simulation  the  fuel  was 
modeled  as  pure  methane  with  r *=  1/4,  whereas  the  value  in  [2]  is  based  on  r = 
1/3.82.  The  latter  value  was  the  result  of  an  actual  analysis  of  the  natural 
gas  used  in  the  experiment. 

The  simulation  was  carried  out  twice,  first  using  the  real  CH^  combustion 
model  and  then  using  the  complete  stoichiometric  combustion  model.  Results  of 
the  simulations  are  presented  in  Figures  28-30  where  the  predicted  mass 
fractions  Cj.  u(^)  plotted  for  k -♦  CH^  , O2  , CO2  , H2O,  CO,  and  H2  . 

As  can  be  seen  from  the  Cj.  ^ plots,  the  GGERM  predicts  that  steady  state 
conditions  are  approached  between  900  s and  1800  s into  the  experiment.  This 
is  consistent  with  the  experimental  procedure  which  involved  steady- state  gas 
sampling  after  approximately  1800  s [2]. 

The  results  for  ^ , Cq2,u«  ^H2,u  (^^®  bottom  plots  of  Figures 

28  and  the  bottom  plot  of  Figure  30,  respectively)  offer  no  surprises.  The 
GGERM  predicts  that  each  of  these  mass  fractions  changes  monotonically  with 
time  from  its  initial  value  to  its  steady  state  value.  As  required  by  the 
GGERM,  note  that  the  latter  steady  state  values  are  consistent  with  the 
results  of  Figures  4,  5,  and  9 or  Figures  16,  17,  and  21  for  <f>  ^ =2.33. 
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The  plots  for  c^q2,u’  ^H20  U’  ^co  u top  and  bottom  plots  of  Figures 

29  and  the  top  plot  of  Figure  30,  respectively)  predict  the  unusual  result 
that  each  of  these  mass  fractions  will  reach  a local  maximum  during  the 
transient  portion  of  the  experiment.  (It  is  noteworthy  that  for  each  of  these 
products  of  combustion  the  predicted  maximum  is  actually  larger  than  the 
maximum  of  the  steady  state  values  for  range  0 < (/)(ss)  < 2.33.) 

The  local -maximum- Cj^  u result  occurs  because  of  the  fact  that  during  the 
transient  portion  of  the  experiment  the  upper  layer  environment  participates 
in  the  combustion  process  by  way  of  the  feedback  flow  at  times  when  the  upper 
layer  already  contains  relatively  significant  mass  fractions  of  the  combustion 
products.  Thus,  the  products  produced  in  the  combustion  process  are  added  to 
the  products  that  already  exist  in  the  layer.  Instantaneously,  the  product 
concentrations  of  the  reactor  output  at  those  times  is  analogous  to  the 
concentrations  of  the  steady  state  upper-layer  environment  in  a reference- [ 2 ] - 
type  of  experiment  where  the  laboratory  atmosphere,  entrained  into  the  fire 
plume  from  the  lower  layer  and  injected  directly  into  the  upper  layer,  is 
contaminated  with  a non-negligible  amount  of  the  products  in  question.  In 
such  experiments,  the  steady- state  upper- layer  mass  fractions  of  these 
products,  as  a function  of  would  exceed  the  steady  state  mass  fractions 

reported  in  [2] . 

While  the  overall  results  of  the  above  GGERM  simulation  are  qualitatively 
plausible,  transient  data  is  not  now  available  to  validate  them  quantitative- 
ly- 


The  Combustion  of  Hexanes  in  Figure- 31  Test  Configurations  - One  of  Two  Hexane 
Experiments  Reported  in  T 3 1 

Consider  the  tests  of  [3].  The  experimental  configuration  for  these  are 
depicted  in  Figure  31.  This  is  similar  to  the  configuration  of  Figure  3a 
except  for  the  fact  that  the  collector  is  a fully  enclosed  space  which 
contains  the  fuel  source.  Air  supplied  to  the  lower  layer  is  brought  in, 
below  the  elevation  of  the  fuel  source,  through  a duct  from  the  outside. 

The  time -dependent  rate  of  air  flow  through  the  duct  is  measured  during  the 
course  of  an  experiment.  This  air  flow  is  not  controlled.  Rather,  it  occurs 
naturally  as  a result  of  cross -duct  pressure  differentials  generated  by 
buoyancy  forces  and  buoyancy -driven  flows  within  and  through  the  collector 
enclosure . 

The  time -dependent  pyrolysis  rate  of  the  fuel  is  measured  by  a fuel -bearing 
load  cell. 

A thermocouple  tree  measures  temperatures  in  the  enclosure  as  a function  of 
time  and  elevation. 

The  GGERM  was  used  to  simulate  the  larger  of  the  two  hexane -burn  experiments 
reported  in  [3].  For  this  experiment  the  thermocouple  data  indicate  that  the 
upper  layer  gases  fill  most  of  the  enclosed  space  relatively  early  in  the 
experiment.  Therefore,  throughout  the  simulation  the  upper- layer  volume  is 
assumed  to  be  fixed  at  1.62  m^ , approximately  ninety  percent  of  the  volume  of 
the  enclosure  above  the  elevation  of  the  fuel  surface.  The  temperature 
history  of  the  upper  layer  is  taken  to  be  the  temperature  measured  at  the 
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elevation  of  thermocouple  4 of  the  thermocouple  tree,  which  is  approximately 
at  the  middle  elevation  of  the  assumed  upper  layer  thickness.  This  tempera- 
ture is  assumed  to  be  a representative  value  of  the  average  layer  temperature. 
The  data  from  thermocouple  4 are  presented  in  Figure  32.  According  to  [10], 
these  data  are  only  valid  up  to  approximately  700  C.  Other  measurements 
indicated  final  peak  and  relatively  steady  temperatures  of  850  C.  In  the 
simulation,  the  upper  layer  temperature  was  calculated  according  to  the  two- 
line-  segment  approximation  of  the  data,  which  is  also  plotted  in  the  figure. 

Measurements  from  the  fuel  load  cell  show  the  time -dependent  hexane  pyrolysis 
rate,  ’ plotted  at  the  top  of  Figure  33.  In  the  simulation,  the  value  of 

this,  which  is  taken  to  be  identical  to  PpygL  flow in  analysis,  was 

calculated  according  to  a three -line -segment  approximation  to  this  data.  This 
is  also  plotted  in  the  figure. 

Data  from  measurements  of  the  duct  air  flow  rate,  plotted  at  the 

bottom  of  Figure  33.  This  air  flow  rate  is  used  to  approximate  flowin’ 

the  instantaneous  air  flow  rate  introduced  into  the  upper  layer  by  entrainment 
from  the  lower  layer  to  the  plume.  In  the  simulation,  the  value  of 
^AiR  FLOWiN  calculated  according  to  a two -line -segment  approximation  to 

the  data.  This  is  also  plotted  in  the  Figure  33. 

Note  that  the  and  approximations  follow  the  corresponding  data  which 

indicate  a step  change  in  these  value  at  t ^ 0 . 

Using  Eqs . (19),  the  time -dependent  value  of  was  calculated  using  the 

above-mentioned  approximations  for  the  time -dependent  values  of  PpuEL.FLOWiN 
and  P^jR  FLOWIN'  is  plotted  in  Figure  34. 

The  GGERM  simulation  of  this  experiment  was  carried  out  twice,  first  using  the 
real  CgH^^  combustion  model  and  then  using  the  complete  stoichiometric  combus- 
tion model.  Only  the  real -combustion-model  simulation  results  will  be 
presented  here. 

Since  experimental  data  was  acquired  for  the  volume  fraction  concentration  of 
O2  CO2 , and  CO,  all  GGERM- computed  results  for  the  product  mass  fractions 
have  been  transformed  into  volume -fraction  concentrations.  The  transformation 
was  carried  out  with  the  assumption  that  the  upper  layer  constituents  include 
only:  the  six  products  of  combustion,  with  mass  fractions  as  computed  from  the 
GGERM;  argon,  with  a specified  mass  fraction  of  0.009;  and  nitrogen.  The 
upper- layer  volume  or  mole -fraction  concentration  of  a product  k will  be 
designated  here  as  Cj^  ^ [vol]  . 

Results  of  the  simulations  for  the  real  combustion  model  are  presented  in 
Figures  35-37  where  the  predicted  volume  fractions  Cj^  u(t)[vol]  are  plotted 
for  k CgH^^,  O2  , CO2  , H2O,  CO,  and  H2  . Included  in’ the  plots  for  O2  CO2  , 
and  CO  are  the  data  for  the  experimentally  measured  volume  fractions  of  these 
products.  Also  included  in  Figures  35-37  are  the  time -dependent  values  of  the 

^kfilACTOUT • 

In  Figure  34  it  can  be  seen  that  the  simulated  step  change  in  fiigEx  *air 

leads  to  ~ o.7  at  t ■=  0.  This  results  in  initial  values  of 

which  are  different  than  the  initial  values  of  the  u ^ ^ simulation, 

the  latter  being  initialized  at  the  concentrations  of  the  ambient  atmosphere. 
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To  the  extent  that  the  [vol ] are  good  approximations  to  the 

in  the  plots  presented  in  Figures  35-37,  the  GGERM  predicts  that  the 
upper  layer  environment  is  basically  in  a condition  of  quasi -steady  state.  As 
can  be  seen  from  the  figures , in  this  particular  experiment  the  GGERM 
simulation  suggests  that  the  quasi-steady  state  approximation  is  reasonable. 
Nonetheless,  for  those  products  of  combustion  where  experimental  data  is 
available  for  comparison,  i.e.,  for  O2  CO2 , and  CO,  the  transient  GGERM 
predictions  are  seen  to  provide  a superior  estimate  of  the  data  than  do  the 
simple  quasi -steady  state  predictions.  To  this  extent,  the  favorable 
comparison  between  the  data  and  the  simulation  results  provides  some  valida- 
tion of  the  GGERM. 

The  following  additional  observations  and  comments  are  offered: 

1.  The  construction  of  the  fj^  functions  for  the  combustion  of  hexane  are 
based  on  the  steady-state  measurements  of  [6].  As  can  be  seen  in 
Figures  10-15,  the  maximum  of  the  4>' s in  those  experiments  was  ap- 
proximately 1.8.  Therefore,  one's  confidence  in  the  present  GGERM 
simulation  beyond  ^'s  of  1.6,  i.e.,  beyond  approximately  t = 150  s, 
should  be  coupled  to  one's  confidence  in  the  extrapolation  of  these  f ' s 
beyond  1.6.  In  other  words,  acquisition  and  use  of  new  data  in  the 
high-0  range  would  lead  to  more  confidence  in  high-4>  GGERM  simulations. 
This  is  of  course  true  for  CH^  and  other  fuels  as  well  as  for  hexane. 

2.  From  approximately  150  s to  the  end  of  the  test,  the  measured  value  of 

Cqz  u reported  to  be  0.003  [10]  whereas  the  predicted  value  is 

approximately  0.03  (see  Figure  35).  The  measured  value  is  not  consis- 
tent with  the  Figure  34  and  the  Cq2  of  Figure  23.  For  example, 
according  to  Figure  34,  at  150  s the  value  of  ~ 2.  According  to 

Figure  23,  this  corresponds  to  ^ = 0.04.  Thus,  for  the  same  4) 

the  Cq2  u results  of  the  steady-state  experiments  of  [6],  on  which 
Figure  23  and  the  present  hexane  combustion  is  based,  is  not  consistent 
with  the  present  quasi-steady  experimental  result  of  [3].  It  has  been 
determined  that  the  discrepancy  is  the  result  of  a non-negligible  depen- 
dence of  fQ2  on  temperature.  In  particular,  it  is  likely  that  for  t > 
150  s the  temperature  of  the  upper  layer  in  the  presently  simulated 
experiment  is  several  hundreds  degrees  C higher  than  temperatures  in  the 
reference  [6]  experiments. 

The  effect  of  similar  levels  of  temperature  difference  on  the  steady- 
state  Cq2’u  CH^ -burn  experiments  can  be  seen  in  Figure  5.  In  this 
figure,  the  data  of  Toner  [5]  were  acquired  at  significantly  higher 
upper -layer  temperatures,  as  much  as  200° C higher,  than  were  those  of 
all  of  the  other  data.  For  ^ > 1,  the  effect  on  Cq2  u significant 
and  completely  consistent  with  the  effects  on  Cq2  u presently  observed 
during  the  burning  of  hexanes.  Presumably,  if  the  hexane  burn  data  of 
[6]  were  acquired  in  a well- insulated  collector  with  upper- layer 
temperatures  in  the  700-850  C range,  the  measured  Cq2  u at  <^  > 1 would 
be  lowered  from  the  order  of  a few  percent  to  the  order  of  a few  tenths 
of  a percent. 

The  expected  qualitative  effect  of  high  temperature  on  Cj,q2,u  ^co,u 

in  the  burning  of  hexanes  can  also  be  obtained  from  the  corresponding 
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high  temperature  effects  in  the  burning  of  CH^ . Thus,  in  Figure  6 for 
</>>!,  the  relatively  high  temperature  of  the  upper  layer  in  the  experi- 
ments of  Toner  [5]  is  seen  to  have  resulted  in  u ten  to 

twenty  percent  higher  than  the  c^q2  u levels  of  the’ relatively  low- 
temperature  experiments  of  [2].  A similar  result  is  noted  for  c^q  ^ in 
Figure  8.  Thus,  for  all  three  of  the  measured  product  concentrations, 
the  high  temperature  effect  in  the  hexane  combustion  model  explains  the 
late-time  discrepancies  between  the  present  GGERM  simulation  and  the 
experimental  data. 

It  is  noted  that  for  an  arbitrary  fuel  of  interest,  any  significant 
dependence  of  the  fj^  on  temperature,  that  may  be  determined  in  steady- 
state  tests,  could  be  easily  incorporated  into  the  GGERM. 


The  Combustion  of  Hexanes  in  Figure -31  Test  Configurations  - A Hypothetical 
Reference- r 31 -Type  Hexane  Experiment 

It  is  important  to  validate  the  GGERM  with  data  acquired  in  experiments  which 
clearly  display  transient  behavior  in  the  development  of  the  upper  layer 
environment,  especially  at  times  of  relatively  large  values  of  c^q . In 
general,  the  experiments  reported  in  [2]  would  be  useful  in  this  regard.  The 
first  experiment  simulated  above  is  an  excellent  case  in  point.  The  refer- 
ence- [2]  experimental  program  is  continuing  [11].  When  transient  data  from 
these  experiments  are  available  it  is  recommended  that  they  be  compared  with 
corresponding  GGERM  simulations. 

It  would  appear  that  experimental  conditions  could  be  chosen  so  that  refer- 
ence- [ 3 ] -t3rpes  of  experiments  would  also  clearly  display  the  transient  upper- 
layer  environment  of  interest.  As  discussed  earlier,  the  test  run  simulated 
above  was  somewhat  limited  in  this  regard  in  that  it  displayed  nearly- quasi - 
steady  behavior  during  the  interesting  times  of  relatively  high  c^q [vol] . One 
reason  for  the  nearly- quasi -steady  behavior  was  that  the  magnitude  of 
^FUEL  FLOwiN  ^®sulted  in  relatively  large  induced  values  of  Pair.flowin* 
the  bottom  of  Figure  33,  the  latter  flow  rate  is  seen  to  have  a characteristic 
value  on  the  order  of  0.07  kg/s.  For  a characteristic  air  density  of  1 kg/m^ , 
and  for  the  assumed  1.62  m^  upper  layer  volume,  this  corresponds  to  a 
relatively  small  characteristic  residence  time,  on  the  order  of  20  s,  for  flow 
exchange  of  the  upper  layer.  (In  contrast  to  this,  the  simulated  reference- 
[2]  experiment  with  a upper  layer  volume  of  3.89  m^  and  Pair.flowin 
0.010  kg/s  has  a characteristic  residence  time  of  400  s.) 

For  the  reference- [ 3 ] test  configuration  it  is  expected  that  a significantly 
reduced  fuel  flow  rate,  leading  to  a corresponding  reduced  air  flow  rate, 
would  lead  to  a display  of  the  transient  behavior  of  current  interest  for 
upper  layer  development.  This  idea  was  investigated  by  using  the  GGERM  to 
simulate  a hypothetical  hexane -burn  experiment  involving  values  of 
Pfuel.flomin  Pair.flowis  tenth  as  large  as  the  Figure  33  values  used 

previously.  The  actual  approximation  for  these  functions  are  plotted  in 
Figure  38.  All  other  parameters  defining  the  experiment  were  assumed  to  be 
unchanged  from  the  previous  simulation.  The  hypothetical  experiment  was 
assumed  to  last  500  s,  where  all  input  parameters  were  assumed  to  be  constant 
for  the  last  250  s.  The  results  of  the  new  simulation  are  plotted  in  Figures 
39-41  in  a manner  completely  analogous  to  the  Figure  35-37  results  of  the 
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previous  calculation.  Favorable  comparisons  between  Cj^  ^ data  from  an  actual 
experiment  (similar  to  this  hypothetical  experiment)  and  Cj^  ^ results  of  its 
GGERM  simulation  would  be  an  important  validation  of  GGERM  capabilities. 


SUMMARY  AND  CONCLUSIONS 

The  Generalized  Global  Equivalence  Ratio  Model  (GGERM)  was  developed  to 
predict  the  generation  rates  of  oxygen,  fuel,  and  other  products  of  combustion 
in  rooms  containing  fires  [1].  The  GGERM  extends  to  general  transient 
conditions  the  global  equivalence  ratio  model  established  during  times  of 
steady-state  in  experimental  studies  involving  two-layer  compartment  fires. 

The  present  work  used  the  GGERM  to  predict  upper  layer  mass  fractions  of 
products  of  combustion  (fuel,  oxygen,  CO,  and  others)  in  these  two- layer  fire 
experiments,  but  during  times  of  transient  response. 

Implementation  of  the  GGERM  requires  knowledge  of  the  fj^  (<^)  function  for  each 
product  k of  combustion  of  a fuel  of  interest.  A fj^  function  predicts  the 
steady- state  generation  rate  of  product  k when  the  fuel  burns  in  Figure- 3- 
type  configurations.  The  fj^  would  be  constructed  from  data  acquired  in 
experiments  of  the  type  described  in  references  [2] -[9].  Generic  rules  were 
developed  for  constructing  the  fj^  from  the  steady- state  data.  The  rules  were 
then  applied  to  data  of  [2]  and  [6].  This  resulted  in  prototype  models  (i.e., 
the  f ) for  the  steady-state  combustion  of  methane  (CH^  ) and  hexanes  (CgH^^^) 
in  Figure- 3 - type  fire  scenarios.  The  products  of  combustion  considered  in 
these  latter  models  include  O2 , CO2 , H2O,  CO,  H2 , and  the  fuel  itself. 

With  the  characterization  of  methane  and  hexane  combustion  in  hand,  the  GGERM 
was  used  to  simulate  transient  upper -layer  environments  in  two -layer  experi- 
ments reported  in  [2]  (methane)  and  [3]  (hexanes).  One  of  these  simulations 
produced  the  unusual  result  of  local  maxima  in  the  mass  fractions  of  CO2 , H2O, 
and  CO.  This  was  explained  by  the  action  of  the  GGERM  feedback  mechanism. 

All  predicted  results  of  the  GGERM  simulations  were  found  to  be  plausible  and, 
where  transient  data  were  available  (upper- layer  volume  fraction  o f O2 , CO2  , 
and  CO  for  the  hexanes  experiment  of  [3]),  predicted  and  measured  results  for 
the  upper  layer  concentrations  of  products  of  combustion  compared  favorably. 
The  most  significant  discrepancies  between  experiment  and  prediction  were 
attributed  to  differences  between  the  (relatively  high)  layer  temperatures  of 
[3]  and  the  (relatively  low)  layer  temperatures  in  the  steady- state  fj^  tests 
of  [6].  For  the  most  part  it  appears  that  these  discrepancies  could  be 
removed  with  use  of  a reference -[ 6 ]- type  data  base  having  a more  extensive 
layer  temperature  variation.  Analysis  of  data  presented  in  Chapter  5 of  [2] 
supports  this  contention. 

Since  comparisons  between  simulation  and  experiment  were  limited,  additional 
validation  of  the  GGERM  under  more  varied  fire  conditions  is  required  before 
it  can  be  used  with  confidence  in  two- layer  zone -type  compartment  fire  models. 

A GGERM  simulation  of  a hypothetical  experiment  in  the  reference -[ 3 ] facility 
was  carried  out  in  order  to  identify  a set  of  experimental  conditions  in  that 
facility  that  would  be  expected  to  highlight  the  transient  character  of  the 
developing  upper- layer  environments.  Calculations  suggest  that  a GGERM 
simulation  which  accurately  predicts  the  experimental  data  from  such  an 
experiment  would  be  an  important  test  of  the  GGERM  capability. 
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NOMENCIATURE 


c(  ss ) N = 1 2 

constants,  Eq . (61) 

, U 

mass  fraction  of  combustion  product  k in  the  upper 
layer 

S,u 

volume  fraction  of  combustion  product  k in  the  upper 
layer 

^k.REACTOUT 

mass  fraction  of  combustion  product  k in  the  reactor 
outflow 

^k  , REACTOUT  ] 

volume  fraction  of  combustion  product  k in  the  reactor 
outflow 

^k,REACTOUT,STOICH 

value  of  Cj^  reactout  stoichiometric  combustion 

c(  SS  ) 

'^k, REACTOUT 

value  of  Cj^  reactout  steady  state  conditions 

p ( s s ) r -TTol  1 

k .REACTOUT 

value  of  Cj^  reactout  under  steady  state  condi- 

tions 

c(  SS  ) 

'-k, REACTOUT, STOICH 

value  of  Cj,  STOICH  steady  state 

conditions 

Eq.  (3) 

f ( STOICH  ) 
k 

value  of  fj.  (<f>)  for  stoichiometric  combustion 

mu 

total  mass  in  the  upper  layer 

*AIR 

measured  mass  flow  rate  of  air  into  reference- [ 3 ] 
facility 

*^HEX 

measured  mass  pyrolysis  rate  of  hexane  in  reference - 
[3]  facility 

*F  E E D B K 

mass  flow  rate  to  the  reactor  from  the  bulk  of  the 
upper  layer 

*F  E E D B K , 1 

Eq.  (27) 

*^F  E E D B K , 3 

Eq.  (28) 

*FLOWIN  [^'fLOWOUT  ^ 

total  mass  flow  rate  into  [out  of]  the  upper  layer 

^ ^F  U E L /^O  XY^STOICH 

stoichiometric  fuel- to -oxygen  ratio 

^k.FLOWIN  t^k.FLOWOUT^ 

sum  of  rates  of  flow  of  all  components  of  product  k 
into  [out  of]  the  upper  layer 
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, U 

mass  of  combustion  product  k in  the  lower  [upper] 
layer 

, U 

net  rates  of  combustion  product  k flowing  to  the  upper 
layer 

^k.REACTIN 

net  mass  flow  rate  of  combustion  product  k into  the 
reactor 

r 

^ ^F  U E L /^O  XY^STOICH 

Tu 

temperature  of  the  upper  layer 

"^U.INIT  f^U.FINAL^ 

initial  [final]  temperature  of  the  upper  layer 

t 

time 

^k 

Eq.  (51) 

^k  , 0 

Eq.  (32) 

^k  , ® 

Eq.  (42) 

^co.N  ^ 1>2,3,4 

constants,  Eq.  (69) 

■^C02 , 1 

a constant,  Eq.  (65) 

“^HZO.N  ^ ^ 

constants,  Eq.  (67) 

\ 

Eq.  (17) 

Eq.  (10) 

o 

Eq.  (56) 

4> 

Eq.  (2) 

,ji(SS) 

value  of  <t>  under  steady  state  conditions,  Eq.  (19) 

*^FUEL 

Eq.  (29) 

STOICH  ) 

v’fuel 

value  of  ^puel  ^^der  complete  stoichiometric  burning 

^k.STOICHl’  ^k,STOICH2 

constants  in  representation  of  fj.  , see  Eq.  (73)  and 
below  Eq.  (71) 

^0  2 

Eq.  (30) 

<^(STOICH  ) 

V'02 

value  of  <f>Q2  under  complete  stoichiometric  burning 

mass  generation  rate  of  combustion  product  k 
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Figure  1. 


The  extended  upper  layer  in  a room  of  fire  origin  which  includes 
the  assumed  negligible -volume  fire  plume. 
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accumulator 

(the  bulk  of  the  extended 


Figure  2.  The  GGERM  representation  of  coobustion  and  flow  dynamics  in  the 
extended  upper-layer  of  Figure  1. 
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Figure  3. 


The  experinental  configuration  used  in  [6] , [7] , and  [8]  (a)  and 
in  [2],  [4].  [5],  and  [9]  (b) . 
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Figure  4.  Conbustion  of  sethane  (CH^).  Dat:a  from  [2].  Analytic 

representations  of  ^eactout * conplete  stoichiometric 

combustion,  Eqs.  (51) -(53)  or  Eqs.  (55)  and  (56)  ( );  and 

for  the  curve-fit  of  Eq.  (58)  ( ). 
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Figure  5.  Conbustion  of  methane  (GU^ ) . Dat:a  from  [2].  Analytic 

representations  of  c^I^^eactoot  * complete  stoichiometric 

combustion,  Eq.  (59)  ( );  and  for  the  curve-fit  of  Eq,  (61) 

( ). 
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Figure  6.  Combustion  of  mte thane  (CH^ ) . Data  from  [2].  Analytic 

representations  of  <^§1 ^reactodt  • complete  stoichiometric 

combustion,  Eq.  (62)  ( ) and  for  the  curve-fit  of  Eq.  (64) 

and  (65)  ( ). 
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Figure  7.  Conbustion  of  methane  (CH^ ) . Data  from  [2].  Analytic 

representations  of  c^lo ^reactodi * complete  stoichiometric 

combustion,  Eq.  (63)  ( ) and  for  the  curve-fit  of  Eqs.  (66) 

and  (67)  ( ). 
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Figure  9.  Conbustion  of  methane  (CH^).  Data  from  [2].  Analytic 

representations  of  c^f^iEACiODi  curve  fit  of  Eq.  (70). 
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Figure  10.  Combustion  of  hexanes  (CgH^^).  Data  from  [6].  Analytic 

representations  of  fc6Hi4  * complete  stoichiometric  combustion, 

Eqs.  (51) -(53)  (71)  ( );  and  for  the  curve-fit  of  Eqs. 

(71)-(73)  ( ). 
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Figure  11. 


Coflibustion  of  hexanes  (CgH^^).  Dat:a  from  [6].  Analytic 
representations  of  fQ2 • fot  complete  stoichiometric  combustion, 

Eqs.  (51) -(53)  and  (71)  ( );  and  for  the  curve-fit  of  Eqs. 

(71)-(73)  ( ). 


39 


N 

O 


CN 

ro 


cn 


Sm' 

O 


1(¥ 


Figure  12.  CcHd>ust:ion  of  hexanes  (CgH^^).  Datia  from  [6].  Analytic 

representations  of  fcQ2 - for  complete  stoichiometric  combustion, 

Eqs.  (51) -(53)  and  (71)  ( );  and  for  the  curve-fit  of  Eqs. 

(71)-(73)  ( ). 
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Figure  13.  Coabustion  of  hexanes  (CgH^^).  Data  from  [6].  Analytic 

representations  of  fH2o  * complete  stoichiometric  combustion, 

Eqs.  (51) -(53)  and  (71)  ( — );  and  for  the  curve-fit  of  Eqs. 

(71)-(73)  ( ). 
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Figure  14. 


Coabustlon  of  hexanes  (CgH^^).  Dafa  from  [6].  Analytic 
representation  of  fj.Q  according  to  the  curve  fit  of  Eq.  (74). 
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Figure  15.  CoHbust:ion  of  hexanes  (CgH^^).  Dat:a  from  [6].  Analyt:ic 

representation  of  fg2  according  to  the  curve  fit  of  Eq.  (75). 
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Figure  16 . Conbustion  of  methane  (CH^  ) . TOP  - Analytic  representations  of 

^REACTOOT • complete  stoichiometric  conbustion,  Eqs.  (51)- 

(53)  or  Eqs.  (55)  and  (56)  ( );  and  for  the  curve-fit  of  Eq, 

(58)  ( ) . BOTTOM  - Corresponding  analytic  representations 

of  f^g^  according  to  Eq.  (23). 
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Figure  17.  CoBbusrlon  of  nethane  (CH^ ) : TOP  - Analytic  representations  of 
*^i^iEACiODi  • complete  stoichiometric  combustion,  Eq.  (59) 

( — ^ );  and  for  the  curve  fit  of  Eq.  (61)  ( — ).  BOTTOM 

Corresponding  analytic  representations  of  £q2  according  to  Eq. 
(24) . 


45 


0.16 


0.14 

0.12 


0.1 

^C02,U  008 
0.06 
0.04 
0.02 
0 

0 0.5  1 1.5  2 2.5  3 3.5  4 


3 


2.5 


2 

^C02  1.5 


1 


0.5 


0 

0 0.5  1 1.5  2 2.5  3 3.5  4 

Figure  18 . Coabustion  of  Bethane  (CH^  ) : TOP  > Analytic  representations  of 
^oi^EEACTOUT  conplete  stoichiometric  combustion,  Eq.  (62) 

( — ),  and  for  the  curve  fit  of  Eq.  (64)  and  (65)  ( ). 

BOTTOM  - Corresponding  analytic  representations  of  f^oz  according 
to  Eq.  (25). 
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Figure  19 . Combustion  of  methane  (CH^ ) : TOP  - Analytic  representations  of 
^io^REACiODT  complete  stoichiometric  combustion,  Eq.  (63) 

( ),  and  for  the  curve  fit  of  Eqs.  (66)  and  (67)  ( )- 

BOTTOM  - Corresponding  analytic  representations  of  fH20 
according  to  Eq.  (25). 
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Figure  20 . Conbustion  of  Bethane  (CH^ ) : TOP  - Analytic  representation  of 
^^o^iEACTODi  according  to  Eqs.  (68)  and  (69).  BOTTCMl  - 
Corresponding  analytic  representations  of  f^Q  according  to  Eq. 
(25). 
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Figure  21.  Conbustion  of  methane  (CH^):  TOP  - Analytic  representation  of 
c^ss^eactoot  according  to  Eq.  (70).  BOTTOM  - Corresponding 
analytic  representation  of  f^2  according  to  Eq.  (25). 
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Figure  22.  Conbustlon  of  hexanes  (CgHj^):  TOP  - Analytic  representations  of 
fpgHiA  complete  stoichiometric  combustion,  Eqs.  (51) -(53)  and 

(71)  ( ),  and  of  the  curve-fit  of  Eqs.  (71)-(73)  ( — ). 

BOTTOM  - Corresponding  analytic  representations  of 
, REACTODT  ^^^ording  to  Eq.  (20). 
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Figure  23.  Conbustion  of  hexanes  (CgH^^):  TOP  - Analytic  representations  of 
fQ2  complete  stoichiometric  combustion,  Eqs.  (51) -(53)  and 

(71)  ( ),  and  of  the  curve-fit  of  Eqs.  (71) -(73)  ( ). 

BOTTOM  - Corresponding  analytic  representations  of 
according  to  Eq.  (21). 
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Figure  24.  CoBbustion  of  hexanes  (CgH^^):  TOP  - Analytic  representations  of 
f(.Q2  coBplete  stoichiometric  comhustion,  Eqs.  (51) -(53)  and 

(71)  ( ),  and  of  the  curve-fit  of  Eqs.  (71) -(73)  ( — • ). 

BOTIOH  - Corresponding  analytic  representations  of 
^olixEACTonT  according  to  Eq.  (22). 
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Figure  25.  CoBbustion  of  hexanes  (CgH^^):  TOP  - Analytic  representations  of 
fgzo  complete  stoichiometric  combustion,  Eqs.  (51) -(53)  and 

(71)  ( ),  and  of  the  curve-fit  of  Eqs.  (71) -(73)  ( — — ). 

BOTTOM  - Corresponding  analytic  representations  of 
‘=AIo1keactoot  according  to  Eq.  (22). 
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Figure  26.  Conbustion  of  hexanes  (CgH^^):  TOP  - Analytic  representations  of 
fj.Q  for  the  curve-fit  of  Eq.  (74).  BOTTOM  - Corresponding 
analytic  representation  of  cJo^izACiooi  according  to  Eq.  (22). 
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Figure  27.  Conbustion  of  hexanes  (CgH^^):  TOP  - Analytic  representations  of 
fg2  fot  the  curve-fit  of  Eq.  (75).  BOTTOM  - Corresponding 
analytic  representation  of  c^f^iKACTOUi  according  to  Eq.  (22). 
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Figure  28.  Result:s  of  GGEEUf  simulation  of  test  1 of  Appendix  B of  [2]  using 
the  real  conbustion  model  ( ) and  the  coBq)lete 

stoichiometric  conhustion  model  ( ).  TOP:  Plots  of 

BOTTOK:  Plots  of  Cg2^||* 
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Figure  29.  Results  of  GGERM  suiulation  of  test  1 of  Appendix  B of  [2]  using 

the  real  conbustion  model  ( — ) and  the  complete 

stoichiometric  combustion  model  ( — ).  TOP:  Plots  of  Ccq2,ji 

BOTTOM:  Plots  of  c^20  o * 
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Figure  30.  Resultis  of  GGERM  sinulation  of  test  1 of  Appendix  B of  [2]  using 
the  real  conbustion  model  ( "■■■“■  ) and  the  complete 

stoichiometric  combustion  model  ( ).  TOP:  Plots  of  c^-q 

BOTTOM:  Plots  of  0^2  u • 
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Figure  31.  Test  configuration  for  the  experinents  of  [2]. 
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Figure  32.  Large -bum  hexane  experiment  of  [3].  Temperature  data  from  the 
thermocouple  located  approximately  at  the  mid- elevation  of  the 
upper  layer  data.  Curve  fit  of  the  data  used  in  the  analysis. 
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Figure  33.  L^rge-buni  hexane  experinent  of  [3].  TOP  - Hexane  pyrolysis  rate 
data  and  curve  fit  of  the  data  used  in  the  analysis;  BOTTOM  - 
airflow  rate  data  and  curve  fit  of  the  data  used 
in  the  analysis  to  approxinate  Ppu^i,  flowii* 
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Figure  34.  Plot  of  foj-  ^he  large-bum  hexane  experiment  of  [3] 
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Figure  35.  Results  of  slanilating  the  large-bum  hexane  experiment  of  [3] 

using  the  real  conbustion  model:  The  GGERM  simulation  ( ) 

and  the  quasi-steatty  approximation  ( ).  TOP:  Plots  of 

bottom:  Plots  of  Cq2  |,[vo1]. 
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Figure  36. 


Result:s  of  sisulating  the  large -bum  hexane  experinent  of  [3] 
using  the  real  combustion  model:  The  GGERN  simulation  ( ■“  ) 

and  the  quasi-stea^  approximation  ( ).  TOP:  Plots  of 

BOTTOM:  Plots  of  0^20. 
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Figure  37.  Resultis  of  sunilating  the  large -bum  hexane  experinent  of  [3] 

using  the  real  conbustion  Hodel:  The  GGERM  sinulation  ( ■■■■  ) 

and  the  quasi-stea<ty  approxination  ( — ).  TOP:  Plots  of 

BOTTOK:  Plots  of 
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Figure  38.  Hypothetical  reference- [3] -type  hexane  experiment.  TOP  - Assumed 
hexane  pyrolysis  rate  used  in  the  analysis;  BOTTOM  - Assumed  air 
flow  rate  used  in  the  analysis. 
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Figure  39.  Results  of  sijBulating  the  hypothetical  reference- [3] -type  hexane 
experiment  using  the  real  combustion  model:  The  GGERM  simulation 

( ) and  the  quasi-steady  approximation  ( )-  TOP: 

Plots  of  u [vol] . BOTT(M:  Plots  of  Cq2  o * 
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Figure  40.  Resultis  of  sinulating  the  hypothetical  reference- [3] -type  hexane 
experiment  using  the  real  combustion  model:  The  GGERM  simulation 

( ■“■■■"  ) and  the  quasi-steady  approximation  ( ).  TOP: 

Plots  of  Cj,Q2  ij[vol].  BOTTOM:  Plots  of  c^20  u * 
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Figure  41.  Results  of  sinulating  the  hypothetical  ref erence - [ 3 ] - type  hexane 
experiment  using  the  real  combustion  model:  The  GGEEIN  simulation 

( ) and  the  quasi-steady  approximation  ( ).  TOP: 

Plots  of  Cj,Q  |,[vol].  BOTTOM:  Plots  of  c^2 
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